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Abstract: We consider an optimal control problem for a vibrating string with multiplicative 
distributive control. Multiplicative control is the coefficient u(x,t) of velocity term y t . Using 
multiplicative control, we bring the state solution to the desired profile under a quadratic cost 
of control. We prove the existence of optimal control. Criterion for characterization of the unique 
optimal control is provided in the form of optimality system. 
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1. INTRODUCTION 

Various engineering applications require prediction of out- 
puts governed by partial differential equations. Many real- 
life applications of optimal control problems with con- 
straints in the form of partial differential equations (PDEs) 
are modeled using hyperbolic equations. Hyperbolic equa- 
tions describe transport processes. Because of their nature 
hyperbolic equations are able to transport discontinuities 
of initial or boundary conditions into the domain on 
which the solution lives or to develop discontinuities in 
the presence of smooth data, these problems constitute 
a severe challenge for both theory and numerics of PDE 
constrained optimization. 

In structural dynamics, dynamic instability created due to 
parametric resonance is an important factor. For example, 
it is believed that instability caused by parametric reso- 
nance is the reason for the famous Tacoma bridge collapse 
in 1940 (Bcrdichersky et al. (1999)). We can control this 
phenomenon by suitable control of the coefficients. Control 
of coefficients has been studied and applied in structures 
governed by elliptic problems (Lurie (1999)), well devel- 
oped literature is available in this case. However, the effect 
produced by bilinear control for hyperbolic problems is not 
studied rigorously (Lurie (1999)-Lurie (2003)). Optimal 
control of systems governed by hyperbolic equations is 
of special importance for the active control of structural 
systems. The field of structural control has been an active 
research area for a number of years. However, most of the 
studies in this area considered specific structures such as 
wings (Wang (1988)), beams (Pan and Hansen (1993)) 
and plates (Bruch (Jr.)). Even though these studies pro- 
vide solutions for many particular cases, the theoretical 
foundations of the subject aimed specifically at problems 
arising in structural mechanics have not received much at- 
tention. Theoretical studies such as the ones in Basile and 



Mininni (1990), Fattorini (1993) considered the optimal 
control problems in abstract settings leaving a gap be- 
tween the theory and applications. In particular, optimal 
control studies relating the theory directly to the solution 
method have been scarce. However, one such application 
is given in Bruch (Jr.) (1995) for a structural vibration 
problem governed by a single hyperbolic equation. An- 
other application is given herein for a vibration problem 
governed by a system of hyperbolic equations. The maxi- 
mum principle (Bruch (Jr.) (2006)) was used to construct 
explicit solutions for an optimal control problem involving 
a distributed parameter structure governed by a system of 
hyperbolic differential equations. Numerical methods for 
solving such problems have not received much attention. 
Sonawane et al. (2013) considered the numerical treatment 
of the optimal control problems for bilinear hyperbolic 
systems. 

In this paper, we study the optimal control problem for a 
dynamic system governed by a hyperbolic equation with 
multiplicative control u(x, t) which is a coefficient of veloc- 
ity term y t . The control u can be interpret as the gain of 
the viscous (motion-activated) damping acting upon the 
string at point x and time t. Bradley, Lenhart and other 
treated the optimal control problem with multiplicative 
control for the plate equation without velocity term in 
Lenhart and Bradley (1994) and with velocity term in 
Bradley et al. (1999). Liang (1999), Lenhart et al. (1999) 
treated optimal control problem with multiplicative con- 
trol without velocity term for wave equation, where as 
Lenhart and Liang (2000) treated similar problem with 
velocity term for wave equation. Belinskiy and Lenhart 
(2001) treated the optimal control problem with multi- 
plicative controls for a wave equation with viscous damp- 
ing. In Lenhart and Liang (2000), a control function is a 
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function of time variable t only, where as in Bclinskiy and 
Lcnhart (2001) it is a function of a space variable x only. 

Let = (0, 1) and T > 0. Consider the following initial- 
boundary value problem for a hyperbolic equation: 

w u = w xx - uw t +vw + /, in Q = il x (0, T), 

w = 0 x = 0, i G (0,T), 

to = 0 x = l, t€ (0,T), (1) 

w = wo(i), t = 0, x G fi, 

tot = i = 0, x G fi, 

where w 0 G H%(Sl), Wl G L 2 (ft), / G L 2 (Q), v G L°°(Q). 
it G Hq(Q) is a multiplicative (also known as bilinear) 
distributive control. We assume that the control functions 
belong to the set of admissible controls: 

U = {ueH 0 l (Q)\0<\u\<M}, (2) 

where M is a constant. 

Consider the objective function: 

J(u) = % [ (w-z) 2 dxdt+^- [ [(u x ) 2 + (u t ) 2 ] dxdt (3) 

1 Jq 1 Jq 

where z G L 2 (Q) is given target function. This objective 
function combines the energy of control and the measure 
of the "closedness" to the desired profile z. We want to 
minimize the function over the set of admissible controls 
U, i.e. to characterize u* such that 

J(u*) = min J(u). 

We shall state and prove the existence of optimal con- 
trol functions. We, then, shall give a characterization for 
control u(x, t) which is a coefficient of a velocity term for 
given v(x, t). The rest of the paper is organized as follows. 
In section 2, we prove the well-posedness of our problem 
and obtain priori estimates. Then, we show that there 
exist optimal control function which minimizes the cost 
functional. In section 3, we derive the optimality system by 
differentiating the cost functional with respect to control 
function u and give the characterization of an optimal 
control through an elliptic variational inequality. 

2. EXISTENCE OF OPTIMAL CONTROLS 

The weak solution of equation (1) is defined as follows: 

Definition 2.1. Given u G U, we say that a function 
w G C([0,T]; Hq(Q)), with w t G C([0, T];L 2 (Ct)), w tt G 
C([0, T]; if is a weak solution of the problem (1) if 

• for any 0 G H%(Q), and a.c. 0 < t < T, 

(w t t,4>) +B[w, 4>\ + / uw t (t>dx= I f4>dx, 
Jn Jn 

• w(x, 0) = wo(x), w t (x,0) = wi(x); 

where (•, •) denotes the duality pairing of £f _1 (£l), 
Hq(SY), and B[w,<f>;t] — f Q w x (j) x dx — J n vwcf)dx is a 
bilinear form. 

Set the notational conventions as follows: 
H = £#(fi) x L 2 (n), 
w = (w,w t ). 

First, we state and prove the existence, uniqueness, regu- 
larity and priori estimates for weak solutions of equation 
(1). We, here assume that the control functions are given 
to study the dynamics of the solutions. 



Theorem 2.1. Assume that w 0 G H 2 (Vt), w\ G Hq(CI), 
u G U n H X {Q), v G C^Q), and /, f t G L 2 (Q). Then 
the problem (1) has a unique weak solution w satisfying 
w e C([0,T};H 2 (n)), w t G C([0, T]; H 1 ^)) and w tt G 
C([0, T]; L 2 (il)). Moreover w satisfies the bound estimate 



sup 

0<t<T 



(\\w{t)\\Hl(n) + \\w t (t)\\ L 2 {n) ) 

+ IkttMllz^acT];//- 1 ^)) 

< C (||/||z,2 ([0;T] . L 2 (Q)) + \\w 0 \\ H i { n) + lki||^(n)) , (4) 

where C is a constant which depends only on Q and T. 

Proof. We will prove the theorem in three steps. In first 
step we will prove the existence of a unique weak solution. 
Regularity of the weak solution will be proved in second 
step 2. And, estimate (4) will be obtained in third step. 

Step 1. We write the equation (1) in the following 
semigroup form 



0 I 

8 2 



dx 2 



0 



w 

m 



o 

-uw t + vw + f 



(5) 



w(0) — w a — 



w 0 

Wl 



Define the operator A : [H 2 (fl) n (fl)} x H%(Sl) -> U by 
Aw 



II I 

d 2 



dx 2 



0 



w 



The domain of A, D(A) = [H 2 {n) n H&(Q)] x H^(Q) is 
clearly dense in W. Equation (5) can be written as 



d_ 
dt 



w(t) = Aw(t) + B(w)(t), 



(6) 



where 



B(w) 



0 

-uw t + vw + f 



Note that A is skew adjoint (Pazy (1983)), i.e., A* = -A. 
Thus D{A) = D(A*), and hence it generates a unitary 
group on H. 

Due to the semigroup form (6), a weak solution should be 
of the form 



w(t) = e At w 0 



oMt-r) 



B{w){-,T)dr, 



where e At represents the semigroup generated by A. 
Define a map S as, 

Sw{t) = e At w 0 + f e A ^B(w)(; r)dr. 
Jo 

We will prove that the map S has a unique fixed point 

Sw(t) = w(t) 
in C([0, Tq];%), where T 0 is to be chosen. 
First, we prove that S is bounded: 
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ll^llc([0,T o ];W) < ll eA **o||c([0,T 0 ];«) 

+ sup f '\\e A ^B(w)(;r)\\ n dr 

0<t<T o Jo 



< sup / ( 

0<*<T 0 



/' 

i JO 



A(i-r)r 



+ sup 

0<t<T 0 Jo 



-U(; T)w t (-,T)+v(-,T)w(-,T)]\\ L 2 {n) dT 

\\e A ^\f(;r)\\ LHn) dr 



+ \\e At w Q \\ c{[ ^ To] , n) 



< 



SUP / \\-u(-,T)w t (-,T)+v(-,T)w(-,T)\\ L 2( n) dT 
0<t<T o JO 

+ sup / ||/(-,r)|| L 2 (n) dr + \\wq\\-h, 
o<t<T 0 Jo 

as ||e A *|| = 1. Since ||u||jyi(Q) < M we obtain 

ll^llcdCTo];^) < ll*o||« + (Coi + Co2)T 0 \\w\\c([0,T 0 ];n) 

+ ^o||/||l~(Q), (7) 
where constant Cqi depends only on M. Hence S bounded. 

Next, we show that S is contractive: For any p, q G 
C([0,T o ];H), we have 



(w m )tt(t)<Jjdx + / u(x,t){w m )to-jdx + B[w m ,<Jj;t] 
Jn 

= f f(x,t)a r (11) 

J LJ 

Multiplying (11) by and then summing over j, we 

obtain 

/ (w m )tt(t)(w m )t(t)dx + / u{x,t){{w m ) t ) 2 dx 
Jn Jn 

+ B[w m ,{w m ) t ;t} = / f(x,t)(w m ) t (t)dx. (12) 
Jn 

Rewriting above equality and using Poincare's inequality 
Evans (1998), we have 



ld_ 
2dt 



/ {{w m ) t ) 2 {t)dx + / (w m ) x (w m ) x (t)da 
Jn Jn 

+ / u(x,t)((w m ) t ) 2 dx 
Jn 



= / v(x,t)w m (t)(w m ) t (t)dx + / f(x,t)(w m ) t {t)dx 



\\Sp - Sg||cr([0,T 0 ];70 < {M v + M)T 0 \\p - g|| C ([0,To];W) 



< 



N 



(8) 



w, 



l {t)dx+ f {{w m ) t f{t)dx 
Jn 



where M„ is a bound for v. 



By selecting T 0 < ( M + M ^ , we conclude that S is contrac- 
tive for t < To. 

From (7)-(8), we conclude that S is bounded and contrac- 
tive on C([0, T 0 ]; "H). Thus, by the contraction mapping 
theorem, there exists a unique fixed point for the map S , „. , , . , , r „ „, 

onC([0 Tj'H) where C 0 is a constant independent of m. Gronwall s 



+ 2 / ((w m )tY(t)dx+- / f\x,t)dx 

< C 0 / (w m ) a; (w m ) a; (t)da; + C 0 / {{w m ) t f{t)dx 
Jn Jn 

+ \ [ f(x,t)dx (13) 



Now, we extend above result to a solution on [To, 2To] by 
selecting new initial value as 



inequality Evans (1998) leads to the estimate 
{{{w m ) t ) 2 {t) + \{w m ) x \ 2 )dx 



sup 

0<t<T 



w Ta = w(T 0 ) e H. 

By applying contraction mapping theorem second time, we < r \ f v * nW . , /* i/,„ n f n \|2,j_ , / f 2 r7r/// 
have a unique solution on C([T 0 , 2T 0 ];U). Repeating the S Gl [ J n [{ m,t> ( ' + J n ^ Wm ' x[V >\ + J Q J a 



process finite number of times, we conclude the existence 
of a unique solution to (1) with w G C([0,T];H). 

Step 2. Let us select a sequence of approximations by 
selecting smooth functions <tj = (Tj(x), (j = 1,2, ■ ■ ■ ) such 
that 



and 



{&j}T=i is a basis of Hq(Q) 



{&j}^Li is an orthonormal basis of L 2 (Q), 



where {aj} are eigenvectors of —A (corresponding to the 
eigenvalues Xj). For integer m, write 

m 

Wm {t)^Y, d T^As) (9) 

where d™-(t) satisfies 

df(0) = (woix),^)^, (df)t(O) = { Wl {x),aj) L 2 

0" = 1,2,- •- ,m) (10) 

and 



Thus we need to estimate the bounds for 

f ((w m ) t ) 2 (0)dx+ f \(w m ) x (0)\ 2 dx. 
Jn Jn 

We have 

m 

\\(Wm)t(0)\\ 2 L 2 in) = H^^l'^O^W^WIU 2 ^) 
m 

< |^](wi,crj)L2 (0) | 

< l|wi||L 2 (n); 



and 



( Wm ) x (0) = ^d™(0)(a j ) x (a ; ), 
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/ (w m )x(0)(w m ) x (0)dx = - / w m (w m ) xx (0)dx 
Jn Jn 

= -J a ^Ew^wj (x>r(°)(^)-(^ dx 

m m 



J'=l 



sup 

0<t<T 



< 



= E A ^ W °' CT ^ 2 - Ikollff^n)- 

For the last inequality above refer Evans (1998). 
Thus, we have 

/ ((w m ) t ) 2 (0)dx+ f \(w m ) x (0)\ 2 dx<C 2 (\\ Wl \\ LHn) 
Jn Jn 

+IKH/fi(n))- 

So, we conclude that 

/ {{{w m )tf{t) + \{w m ) x \ 2 )dx 
Jn 

C 3 [\\wi\\ LH n) + IKHffi(n) + \\fh'(a)] • (14) 
From (11), we obtain 

/ (w m ) tt (t)<j)dx + / u(x,t)(w m ) t (f)dx + B[w m ,4>;t] = 
Jn Jn 

/ f{x,t)cj)dx 
Jn 

for any 0 G ff^(fi). Thus, 

+|B[io rn ,0;i]| 

+ / f(x,t)<j)dx . 
Jn 

Poincare's inequality and inequality (14) together gives 



/ (w m ) tt {t)4>dx 


< 


/ u(x,t)(w m ) t <f)dx 


Jn 




Jn 



/ {w m )tt(t)<j>da 
Jn 



for any <j> S iJ^fi). Passing to the limit on a subsequence 
as m — » oo, we obtain u> e £ 2 ([0, T]; i? 1 (fi)) such that 

w m — > w weakly in L 2 ([0,T}; H 1 ^)), 
(w m )t — > w t weakly in L 2 (Q), 
(w m ) tt — ► «;„ weakly in L 2 ([0,T]; (ff 1 )*^)). 

Step 3. Consider the sequences {wo n } C i? 2 (f2), {wi n } C 
i?o(0) and {u„} C {/ such that 

w 0n — > w 0 strongly in Hq(Q,), 
Wi n — ► w i strongly in L 2 (f2), 
u n — > u strongly in L 2 (Q), 

Let w n be the weak solution of (1) corresponding to the 
initial values u>o n , w\ n and control u n , then w n satisfies 
the regularity of previous step. Multiply the equation (1) 
by (w n )t and integrate over Q t — ft x (0, t) with 0 < t < T, 
we obtain 



0 = / [{wn)tt{w n )t + u((w n )t) 2 - {w n ) xx (w n ) t } dxdr 
jQt 

[vw n (w n ) t + f(w n ) t ] dxdr 



J 

Jq, 



Qt 



ld_ 
2dt 



(((w n ) t ) 2 + \(w n ) x \ 2 )+u((w n ) t )' 



J 

JQt 



dxdr 



[vw n (w n ) t + f(w n ) t ] dxdr, (15) 



where for 0 < t < T, w n (Q,t) = w n (l,t) = 0, {w n ) t (0,t) = 
(w n )t(l,t) — 0. From (15), we obtain 

/ [((W n ) t ) 2 + \( Wn ) x \ 2 ]dxdT 

Jnx{t} 

= / [(w ln ) 2 + \(w 0n ) x \ 2 ] dx 
Jn 

+ 2 [-u((w n ) t ) 2 +vw n (w n ) t + f{w n ) t ] dxdr, 
Jqt 

< \\ w ln\\ 2 L2 (n) + \\(W0n)x\\h(n) + \\f\\h(Q t ) 

+ \\{wn)t\\h(Q t ) + C 5 (\\w n \\ 2 L 2 {Qt) + \\(w n )t\\ 2 L 2 (Qt) ) (16) 
where constant C 5 depends only on M. Using Poincare's 
inequality, we obtain 



/ [(K)*) 2 + K) 2 ] dx 

Jnx{t} 

< C 6 ( [{{w n ) t f 
JQt 



+ (w n ) 2 ] dxdr 



+ C 7 (\\w ln \\ 2 LHn) + \\(w 0n ) x \\ 2 L 2 (n) + \\f\\h iQt) ) (17) 

where constants Cq, C-j are independent of w n . Applying 
Gronwall's inequality to (17), we obtain 

/ [(K) t ) 2 + K) 2 ] dx 

Jnx{t} 

<(l + CeTe c ° T )C r (\\w ln \\ 2 LHn) 

+ \\^onU 2 LH n) + \\f\\h ( Q t) ) (18) 
Letting n — > oo and using inequalities (16), (18), we obtain 

+ \\(wUt)\\h ( n) + Mt)\\ i 2 (n) 
< c 8 (|H||i a(n) + W(woUh m + \\f\\h {Qt} ) (19) 

where constant Cg depends on f2 and T. Thus, 
sup (\\w(t)\\ 2 Hl{Q) + \\w t (t)\\ 2 LHn) ) 

0<t<T 0 ' 

< C 8 (|kl||!2 ( n) + \\w 0 \\ 2 H i m + ||/||!2( [0> T];L2(O)))- (20) 

Next, we have 

(w t t,4>){t) = - I [w x (j) x + uw t <j) - vw<j) - f4>](t)dx 
Jn 

for any <j) € H^(n). Let \\(f)\\ H i (n) = 1, then 

\(w tt ,cf>)(t)\ <C 9 (|| Wt ||| 2(0) + ||w|| 2 Hol(0) ) + ||/||| 2(0) ) (21) 
where constant C 9 depends on M. Using (21), we obtain 

ll W «llL2([0,T];H-i(n)) < Cio / (\\wt\\l2( n) + \M 2 H^n))dt 

J 0 

Il/lli2 ( n)^ 

< CuTdl^lH.^ + \\ W f Hlisl) ) + \\.f\\h (n) ) (22) 
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where Cio, Cn are constants which do not depend of u 
and v. 

Inequalities (20) and (23) together imply the inequality 
(4). ■ 

Now, we prove the existence of optimal controls. 

Theorem 2.2. (Existence of optimal controls). There exist 
an optimal controls u* € U which minimizes the objective 
functional J(u) for u e U. 

Proof. Let {u n } e U be a minimizing sequence such that 



lim J(u n ) = inf J(u). 

ra-s-oo ueU 



Let be a solution of equation (1) with control u n . 
Denote w n — w(u n ). By Theorem 2.1 we have 



sup (\\w n (t)\\ H , {n) + \\(w n ) t (t)\\ LHn) ) 
o<t<T v 0 ' 

+ \\(w n )tt(t)\\m[0,T];H-HQ)) 
< C (||/||L2([0,T];L2(n)) + ||«>o||iji(n) + ll w l|U 2 (fi)) , 



where C is a constant. 

On a subsequence, by weak compactness, there exists w* 
in C{[0,T];H^(n)) such that 



w n — > weakly in L°°([0, T]; ^(fi)), 

(w™) t — > ( W *) t weakly in L°°([0, T]; L 2 (f2)), 
{w n ) tt — > («;*)« weakly in L°°([0, T]; _ff _1 (0)), 

u" — ► u* weakly in if 1 (Q), 
(« n ) x — > (u*) x weakly in L 2 (Q), 
(u n ) t ^(u*) t weakly in L 2 (Q). 



Using compactness result from Simon (1987), we have 
w n -> w* strongly in L°°([0,T]; L 2 (Q)). By definition of 
weak solution, we have 



r \{w n ) tt ,4>)dt 

Jo 

[{w n ) x <j> x + u n {{w n ) t )cp - w n 4> - f<f>]dxdt. 



T A 



0 JO 



{u n {w n ) t 4>-u*{w*) t 4>) dxdt 



f [ \(u n -u*)(i>\ \{w n ) t \dxdt 



+ 



< 



0 JO 
T A 



0 JO 
T A 



u*<j)((w n ) t ~ {w*) t )dxdt 

\ 1/2 



// 

JO JO 



u n -u*)(f)\ z dxdtj 

\ 1/2 



' dxdt 



+ 



[ [ u*<P{{w n ) t -{w*) t )dxdt 
Jo Jo 

\ 1/2 



Jo lo 



<Ci 



\u n -u*\ 2 dxdt 



T A 



0 Jo 



1/2 



\(w n ) t \ 2 dxdt 



T A 



0 JO 



u*4>{(w n ) t - {w*) t )dxdt 



Oasn^oo, 



where C\ depends on the L°° bound on <f). 

The first term converges, since u n — > u* weakly in H 1 (Q) 
and {{w n )t} is uniformly bounded in L 2 (Q). Also, u*<f> G 
L 2 (Q) and (w n ) t — ► {w*) t weakly in L°°([0, T]; L 2 (Q)), 
the second term converges. Thus, 

u n (w n ) t — -»■ u*(w*) t weakly in L 2 (Q). 

Passing to the limit as n — » oo in the weak formulation of 
w n , we obtain 

f \(w*) tu <P)dt 
Jo 

= ~ I [{w*) x <P x +u*{{w*)t)(t)-w*<t>-f<P\dxdt. 
JQ 

Thus w* = w(u*) is the solution of state equation (1) with 
controls u* and v*. Since 

J(u*) = ? / {w* -zfdxdt 



+ § / [(K),) 2 + ((^) t ) 2 ]^- 



'Q 

Using lower-semicontinuity of L 2 norm with respect to 
weak convergence, we have 

J(u*) < lim f / [(K),) 2 + (K) t ) 2 ] dzdi 

+ lim % / (w™ - zfdxdt 

n—>-oo 2 Jq 

< lim inf J(u") 
= inf J(u). 

u£!7 



We have, 



Thus, we conclude that u* is an optimal control. 
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3. CHARACTERIZATION OF THE OPTIMAL 
CONTROL 

By differentiating the objective functional J(u) with re- 
spect to the control u, we derive the optimality system. 
As w = w(u) is involved in J(u), we must first prove the 
appropriate differentiability of the mapping 

u — > w(u). 

Lemma 3.1. The mapping 

u e U -> w(u) = (w(u),w t (u)) e L 2 ([0,T];H) 

is differcntiable in the following sense: 

w(u + eh) — w(u) ^ ^ 

e 

weakly in L 2 ([0,T];V.) as e — > 0, for any u satisfying 
u + eh e U for small e and /i e Hq(Q). 

Moreover * = (*,* t ) e i 2 ([0,T];"H) is a weak solution 
of the following problem: 

*t = t ra + u*t + - to t , in Q, 

* = 0, i = 0, t£ (0,T), 

* = 0, x = Z, te (0,T), (23) 

* = 0, i = 0, x e O, 
* t = 0, t = 0, .t e O, 

where w = w(w). 

Proof. Let w e = w(u + eZ) and u> = u>(w), then (w e — w)/e 
is a weak solution of 



w — w 



w — w 



w — w 



+v 



w — w 



hw e t , in Q, 



e 

w e — w 



= 0, a; = 0, t e (0,T), 



e 

w e — w 
e 



0, a; = Z, t e (0,T), 
= 0, f = 0, ir e n, 
= 0, t = o, xen. 



Using Theorem 2.1, we have 



sup 

0<t<T 



(w e — w) 


+ 


wt ~ w t 


e 




e 



L 2 (f2) 



+ 



Wit ~ w t t 



L 2 ([0.T]iH- 1 (n)) 
< C||/lW t C || L 2([ 0 ,T];L 2 (0)) 

< Ci||/i|| L oc (Q) supO < t < T|| Wt e || L2(0) < C 3 

where C3 is a constant. Hence on a subspace, by weak 
compactness, we have 

Wtt ~ WU — ► *« weakly in L°°([0, T]; if" 1 (ft)), 



wf- ir, 



UT — w 



* t weakly in L°°([0, T]; L 2 (fl)), 
> * weakly in L°°([0,T}; H^fl)). 



w tt 



In ( 



,0) + 



n 

w e — w 



w 



4>dx + 



h 



c dx 



w 



I hw e t 4>dx 
J a 



for any <j) <G H£(Q), a.e. 0 < t < T. 

Letting e — > 0 and using w e — > w, w\ — > w t weakly in 
L 2 (fl), we obtain that ^ is the weak solution of (23). ■ 

Now we prove the necessary conditions that characterize 
an optimal control. 

Theorem 3.1. Given an optimal control u and correspond- 
ing w — w(u) = (w,w t ), there exists a weak solution 
P = (PiPt) in H to the adjoint problem: 

Ptt = Pxx + {up)t +vp+ (w - z) in Q, 
p = 0, x = 0, te (0,T), 

p =0, x = l, te (0,T), (24) 

p = 0, t = t, x e n 
Pt = 0, t = T, x e O. 

Furthermore, u satisfies the following variational inequal- 
ity: 

a 

mm(max(-(u tt + u xx + -pw t ),u- M),u + M) = 0, (25) 
where w — w(u). 

Proof. Let u G U be an optimal control and w — w(u) 
be the corresponding optimal solution. Let u + eh G U 
for e > 0 and w e — w(u + eh) be the corresponding 
weak solution of equation (1). We compute the directional 
derivative of the objective functional J(u) with respect to 
u in the direction of h. Since J is supposed to attain its 
minimum for u, we have 

J(u + eh) - J(u) 



0 < 


lim 








lim 






+ 


lim 




e-+0+ 


+ 


lim 








«/ 




JQ 



Q 



By the definition of weak solution, we have 



[(w £ - z) 2 - (to - z) 2 ]dxdt 
+ lim j|- / [(u t +e/i t ) 2 - (u t ) 2 ]da;d4 



where ^ is as defined in Lemma 3.1. We need to prove 
the existence of the solution to the adjoint equation (24). 
It does not follow from the existence and uniqueness of 
solution of the state equation (1) as adjoint equation 
contains the term (up) t . 

We write the adjoint equation in semigroup form to as 

Jt m = Am+ {(up) t + vp + w-z 
where A generates a unitary group on H. 
We expect a solution of (24) of the form 

(up) t + vp + w - z 



P(t) 



0 A(t-s) 



ds. 
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Integration by parts gives, 

i-T 



Mt-s) 



0 

vp + w — z 



ds + 



0 

up{t) ) ■ 



Now consider the map 
Fp(t) 



o \\,.fo 

ds - 



/ , v V vp + w ~ z ) ) ^ V up ^ J ' 

Using technique similar to the proof of Theorem 2.1, step 
1, we can prove that the map F has a unique fixed point 
on [0,T], which guarantees the existence and uniqueness 
of weak solution to the adjoint equation. 

Assume that p is the weak solution of the problem (24), 
i.e., p satisfies 



(Ptt,<f>){t) + / p x <t> x {t)dx= / {up)t(j)dx+ I vp<j)dx 
Jn Jn Jn 

+ (w- z)(j>(t)dx 
Jn 

for any (f> € H^(Q), and a.e. 0 < t < T. Then 

0 < a I ^>(w — z)dxdt + 8 / [u x h x + u t h t ]dxdt 
JQ JQ 

= a \( Ptt , tf)(t)+ / p x * x (i)da;- / (up) t * \ dtdx 
Jo I Jf2 Jn J 

a J {y w?5 *^ a '| ^ ^ y [ u xh x + u t h t ]dxdt 
~ a J ^yj —hwtpdx^ dt + f3 J [u x h x + u t h t ]dxdt 



J 

Jq 



{—ahw t p + f3{u x h x + u t h t ]} dxdt. 



The last expression gives the variational inequality (25). 
■ 

We obtain the optimality system consisting of two hyper- 
bolic PDEs and an variational inequality, (1), (24), and 
(25). 
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